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SUMMARY 


A method has been developed to predict passenger satisfaction with the 
ride environment of a transportation vehicle. This method, a general approach, 
has been applied to a commuter-type aircraft for illustrative pxirposes. Here 
the effect of terrain, altitude and seat location were examined. The method 
predicts the variation in passengers satisfied for any set of flight condi- 
tions. In addition several non-commuter aircraft were analyzed for comparison 
and other uses of the model described. The method proposed has advantages for 
design, evaluation, and operating decisions. 


INTRODUCTION 


The purpose of this paper is to provide a method of assessing passenger 
satisfaction with the ride quality on transportation vehicles. The method 
is applicable to both existing systems as well as future ones, and can be used 
for evaluation, design and decision making. Basically it relates the 
environment in which the vehicle must be used and the performance character- 
istics of the vehicle to determine the probability of satisfying the 
passenger. 

This analysis is based on previous work by the authors in assessing 
vehicle ride quality for the air mode. In refs. 1 and 2, a model of passenger 
ccmfort and satisfaction with a ride as a function of the motion of the 
vehicle was developed. This model coupled with standard techniques for 
analyzing a vehicle's motion allow us to examine such variables as: vehicle 

type, input forcing functions, operating characteristics, etc. 

The method will be applied to commuter-type aircraft and variations in 
passenger satisfaction due to terrain, altitude, equipment and location 
in the vehicle described. Other uses of the technique are also suggested. 
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SYMBOLS 


a acceleration 

C comfort rating 

f Joint probability density function 

L length scale for turbulence spectr\an 

S percent of passengers satisfied with ride 

Vo velocity of vehicle 

turbulence gust velocity in transverse direction 

V turbulence gust velocity in vertical direction 

S 

y mean 

p correlation coefficient 

a standard deviation for accelerations, rms for turbulence quantities 

(P power spectral density 

(0 frequency 

« I =.! ^ang;ular velocity 

_ ;rta=-V'2' i.-r-ir 

Subscripts: 

a transverse rms acceleration 

y 

a vertical rms aicceleration 


z 

X longitudinal direction 

y transverse direction 

z vertical direction 


METHOD 


Description 

The method of analysis is shown in figure 1. A vehicle forcing function is 
converted into motion cues to the passenger using the appropriate transfer 
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functions for the system being analyzed. Typical forcing functions are illus- 
trated in table I along with the important properties of the transfer functions 
for several vehicles. 


Table I 


Vehicle 

Primary Forcing Function 

Characteristics of 
Transfer Function 

Airplane 

Atmospheric Turbulence 

Aerodynamics, Mass 
Properties 

Train 

Rail Profile 

Suspension System, 
Wheel Configuration 
Mass Properties 

Bus, Automobile 

Road Surface 

Suspension System, 
Wheel Configuration 
Mass Properties 

Ship 

Sea Surface 

Hydrodynamics, Mass 
Properties 


In most cases the vehicle engines also contribute to the motion experiences (e.g. 
vibrations) however their amplitudes and frequencies compared with the primary 
forcing function shown above are usually negligible. 

Vehicle functions genereLLly depend on frequency; thus, both aS^litude and “ 
frequency information of the input; or, the input power spectrum is necessary 
for the analysis. In addition the inputs can and usually are statistib'^l^" ' 
varying quantities so that a probability density fimction for each of the in- 
puts is necessary. In fact, as will be seen below, the method ' described sLLlowl 
for isolating components of the forcing function which contribute most to 
passenger dissatisfaction. In some cases this information may be used to 
find ways to improve the ride environment (e.g. treatment of roadways or active 
ride smoothing on the vehicle). jj 

Vehicle motion can take the form of velocities, accelerations and rates of 
change of acceleration in each of six-degrees-of-freedom. Not all of them are 
appropriate for the ensuing analysis and only those needed in the subjective 
transfer function must be determined. In general the passenger's comfort will 
be functionally related to the motion parameters of angular velocity and linear 
acceleration and their derivatives 


C s f(a ,a ,a ,(o ,b> ,(d ,d ,i ,a ,d) ,d> ,d) ) 
X y z -X -y -z’ x y z -x -y -z 


( 1 ) 
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where C is the subjective comfort rating, a ,a ,a linear accelerations in the 

X y z 

longitudinal, lateral and vertical directions respectively, and o)^, 0 )^, 0 )^ angu- 
lar velocities about the longitudinal, lateral and vertical axes respectively. 
The • denotes a time derivative, thus a^ is the longitudinal jerk, etc. The 

comfort model can be a simple function of rms motion variables through a more 
complex frequency-dependent psychophysical model (see e.g. ref. 2). In the 
"real world" other factors also contribute to the passenger's comfort (e.g. 
noise, temperature, etc.) however these will be neglected here. A more complete 
analysis should include them. 


The mathematical procedure for arriving at the comfort rating is straight- 
forward but somewhat tedious to perform. The joint probability density function 
for the motion variables, f(a^,a^,a^, ...) is integrated over motion space to 

arrive at a probability function for the passenger fe comfort level. That is, 
the probability that the comfort rating is less than or equal to some value C 
is given by 


C'-C 

0 ) 

P(C < c) = / "^/ 

d) =0 


C'-C 



a =0 

X 




(I) )da da da . . . d(l) (2) 
-z X y z -z 


where is the value of the associated motion variable given by the comfort 

equationt each one being eliminated as the integration progresses. Since the 
motion can vary with location in the vehicle, the above analysis must be re- 
peated at each station of interest. 

The last step in the analysis relates the derived comfort rating to a value 
judgement. This value judgement is taken to be passenger satisfaction with the 
ride which is related to comfort rating. The percentage of passengers satis- 
fied, S, is a simple function 


s = f(C). (3) 

Thus for any comfort rating the value judgement transfer fxmction transforms C 
to S by the above equation. The actual decision process is much more complex, 
being dependent on other variables as well as competing modes. These have been 
neglected in this analysis, assuming that if a passenger were dissatisfied a 
sufficient niomber of times he wo\ild seek an alternate means of reaching his 
destination. 

The remainder of this paper will apply this method to a particular vehicle 
type — commuter aircraft — however it is important to note that the method is by no 
means restricted to this mode. At the present time this is the only mode for 
which data were available. 
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Application to Connnuter-Type Aircraft 


Input Forcing Function 


For aircraft the input forcing function is atmospheric turbvilence, which 

can he characterized by velocity power spectra in all six-degrees-of- freedom, 

longitudinal, lateral, vertical, pitch, roll, and yaw. However, since previous 

work (refs. 1 and 2) has shown that the comfort models require only vertical 

and lateral linear accelerations, only these components of the turbulence field 

will be considered. The amplitude probability as well as the frequency content 

are functions of terrain, altitude, and weather. Typical examples are shown in 

figure 2, where the variation in vertical, O , and lateral, a , rms gust in- 

w V 

g g 

tensity is seen versus altitude for moimtain terrain (refs. 3 and 1+). Similar 
CTorves are available for water and flat terrain. The power spectra for these 
are given by a Dryden model 


<t>(u) 


L 

ttVo 


1 + 3 (^)^ 


(h) 


where Vo is the aircraft velocity, o), the frequency, O the rms gust intensity 
and L, the length scale which is a function of altitude (ref. U). A typicsQ. 
power spectrum is shown in figure 3. 


Vehicle Transfer Function 


Aircraft transfer functions are a fvinction of aerodynamics and mass proper- 
ties and can be found in many references (see e.g. ref. 5)* Here we assume a 
rigid body model (no structural bending) and neglect gyroscopic effects. The 
particiilar vehicle first considered is the deHavilland Twin Otter aircraft, 
which was selected because of the abundance of data available concerning its 
aerodynamic characteristics. It is regrettable that functions for the air- 
craft suitable for potential use in the commuter market are not readily avail- 
able. 


Motion Spectra of Vehicle 

The outputs of interest for the comfort model to be used below are the rms 
accelerations in the vertical and lateral directions. These can be obtained 
by integrating their power spectral densities over frequency space which are 
given by 


z 

<t> (w) 

a 

y 


a 2 

f- m 


g 


a 2 

4>(w) 

g 


(5a) 

(5b) 
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where <|)^ (w) (u) are the power spectral densities for vertical and lateral 


z y 

accelerations, and 


g 


and 


a 

_z: 

▼ 

g 


are the transfer functions for these accelera- 


tions relating them to the turbiaence field. For the Twin Otter, the rms ac- 
celeration cumulative probability distribution is shown in figure k for a typi- 
cal case, lypical spectra for these accelerations are given in figure 5.* As 
can be seen in figure 6 the acceleration in the vertical direction closely 
approximates a normal (Gaussian) distribution. The same behavior can be seen 
for actual flight data in reference 6. Transverse acceleration behaves similar- 
ly . This allows us to write the probability density functions for each separ- 
ately and for both combined using a normal distribution. From flight data 
(refs. 2 and 6) the cross correlation between vertical and lateral accelerations 
is 0.8 thus the joint probability distribution f\mction is given by 


f(a ,a )-r— 
y z 2Tra 


a 
a a 
y z 


•yiT 


:exp< - 



2p(a -Ma 

y z 

a a 
a a 
y z 


( 6 ) 


^ere y ,y are the mean rms accelerations, a ,a are the standard devia- 
a a. 8> 8. 

y Z y z . 

tions of rms accelerations, and p is the correlation coefficient between acceler- 
ations. The values for the y’s and a's for different terrain, alij^itude and 
vehicle location can be fo\md by computing values of the motion, variables for 
= ’-.^^anci -84 respectively. ^ 

Subjective Transfer Function , 

A subjective comfort model has been developed (ref. 2) based on extensive 
field data taken on commercial airlines (refs. 6, 7, 8) and in-flight simulator 
(ref. 9) experiments. This model relates the subjective comfort response to 
rms vertical and transverse accelerations in g's as 


C = 2 + 11.9a 


+ 7.6a when a ^T.6a , 
y i z 


(7a) 


and 


C = 2 + a + 25a when a < 1.6a , (7b) 

z ' Y z y* ' ' ' 


*The contribution to the rms acceleration for the vertical direction can be di- 
vided into two ^equency regimes-below and above 1 rad/sec. The rdgion above 
1 rad/sec contributes 88 percent of the totcJ. power and is thus more important 
in determining comfort. 
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>^ere C is restricted to values 2 through with the following desciriptors: 


C = 2 comfortahle 
= 3 neutreJ. 

= k uncomfortable 
= 5 very uncomfortable. 


For motions in which vertical acceleration dominates (i.e. 


\ > 1.6ay), 


subjective Judgments lean more heavily toward the vertical stimulus, however 
transverse acceleration is more important otherwise. For pure motion in either 
direction these models predict twice the sensitivity to the transverse direc- 
tion compared with the vertical direction. 


Comfort Determination 


Using equations 7 we compute the comfort rating corresponding to any givcsi 
vertical and transverse accelerations. However the accelerations cure described 
by the Joint probability distribution function (given in eqviation 6). Thus the 
probability of exceeding a given comfort level C is obtained from equation 6 
using equations 7 to describe the integration space as 


a (C'-2-7.6a )/ll.9 

,yo . y 

P(C > C’) = / / f(a ,a )da da 

0 I.L y ^ y * 
y 


a (C»-2-a )/25 
Zq z 

+ / / f(a a )da da . 

0 a /I. 6 y z z y 
z 


( 8 ) 


For the Twin Otter aircraft the first, y, suid second, o, moments describing the 
probability distribution f sure given in table II, for the center of gravity o^ 
the aircraft. Similar data have been generated for other positions withfn^the 
craft . 


Table II 


Terrain 



Mountain 

„ Water 

Flat 1 

Altitude 

152 m 

3,0U8 m 

152 m 

3,OU8 m 

152 m 

3,0U8 m 

(g's) 

.055 

.035 

.019 

.012 

00 

0 

• 

.031 

y^^ (g‘s) 

.015 

.009U 

.0051 

.0032 

.013 

.0082 

a (g’s) 

cL 

-4* 

CVJ 

o 

• 

.015 

.0052 

.0033 

.019 

.012 

z 

a (g’s) 
a 

y 

.0066 

.001^1 

.OOllt 

.00088 

.0051 

.0032 
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Equation 8 is numerically integrated to determine the cumulative probability 
distribution for each case of interest (i.e. terrain, altitude, seat location). 
A typical result is shown in figure 7 which illustrates the variation due to 
terrain for a fixed altitude and location within the aircraft. Thus for this 
case there is a 90^ probability that the subjective comfort rating will be less 
than 2.3 for flight over water, less than 2.75 for flight over flat terrain, 
and less than 2.9 for flight over movuatain terrain. Stated ea.tematively the 
probability that the comfort rating of the ride will be less than or equal to 
2.5 is 100^, ^h%j and h2% for flight over water, flat and mountainous terrain 
respectively. Similar comparisons can be made for any set of conditions or to 
compare different aircraft for a single set of conditions. 

Value Transfer Function 

The calculated comfort judgments must now be related to a more value- 
oriented variable. We choose as this queuatity the percentage of passengers 
satisfied with the ride, that is, the fraction of passengers who would willingly 
take another flight at least without hesitation. This quantity has been deter- 
mined in previous work (ref. l) to be related to the subjective comfort rating 
as shown in figure 8. As can be seen, from a statistical point of view, there 
are approximately ^% of the passengers who will not be satisfied with the ride 
environment even when the ride is rated comfortable by most of the passengers. 
This is seen more clearly when examining distributions of passenger responses 
(see e.g. ref. 2). 


This transfer function, figure 8, has been applied to data on subjective com- 
^ort responses, to obtain the probability of satisfying a given percentage of 
the passengers. Typical graphs are given in figures 9 and 10 for the aircraft 
center of gravity and an extreme aft seat location as a function of terrain and 
altitude. As -an ex^plp of , using these graphs, they ipdicatW_thae''^tHhe' depter 
of gravity there is. a, 1+5^ probability , pf/satisfying at' leaH ti>*'.p4's^3> ‘ 

gers flying at 152 m over mountain terrain,, while there is an 8UjS probability” , 
of satisfying the same number of passengers flying at 3,0U8 m over mountain ' 
terrain. Similarly over the same terrain at 152 and 3,0U8 m respectively 
there is a 36^ and ^Q% probability of satisfying 85^ or more of the passengers 
at an aft seat location. This illustrates that a) the aft seat locations are 
less comfortable than those near the center of gravity, and b) flying at higher 
altitudes increases the' probaljii^ii'^ o^ satl^fVinV'i>ass^^ Tlius^the^more 

conservative approach woiild be to design to the low altitude, aft seat loca- 
tion resxilts. 




^:7J8 ^hLT;tLtiK 


Comparison to Other Aircraft 

Several other aircraft have been analyzed using the method described. 
Transfer functions were obtained from references 10 and 11. The aircraft are 
the Breguet 9^1, Douglas DC-8, Cessna l82, and an externally blown flap (EBP) 
aircraft still in the design stage. These aircraft have the following char- 
acteristics: 
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Aircraft 

Weight (kg) 

Approximate No. ( 

DC-8 

91,000 

200 

Cessna 182 

1,360 

k 

i^eguet 9^1 

20,000 


EBF 

122,000 

270 


Figure 11 illustrates the variation in percent satisfied hy aircraft type 
for cruise at 3»0U8 m altitiide over mountain terrain. As is seen the DC-8 
is the best aircraft emd the EBF the vorst. 


Applications of the Method 

The method desczdbed can be used to assess the satisfaction of passengers 
with the ride environment of a given vehicle. In addition it can be vised to 
perform sensitivity analyses of the effects of vehicle variables, through varia- 
tions in the vehicle transfer function and of input veuriables through the forcing 
function. These csua be used to determine maximum design payoffs in the case of 
the vehicle or operating conditions and surface reqvilrements (for roadways /rail) 
in the case of the forcing fvinction. 

Another application would be to incorporate an optimization routine and 
use the method Inversely to determine the optimum design under engineering con- 
straints for a desired satisfaction level. 

Lastly, the method cam be applied to validation studies of models of comfort 
and/or satisfaction by testing over a wide ramge of conditions with a limited 
set of field data. This would be acconqplished by inserting the appropriate 
tramsfer function to replace those described above. 


Conclusions 

A method has been developed to predict passenger satisfaction with the ride 
envlTOnment of a transportation vehicle. This method, a general approach, has 
been applied to a commute 3 ?-type aircraft for illustrative purposes. The effects 
of terrain, a^.tltude amd seat location were examined. The method predicts the 
variation in passengers satisfied for any set of flight conditions. Several 
non-commuter aircraft were elLso amalyzed for comparison and other uses of the 
model described. The method proposed hem advamteiges for design, eved.uation, and 
operating decisions. 
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Figure 4.- Acceleration probability distribution. 








Figure 6.- Comparison of normal acceleration distribution 
with Gaussian distribution. 




Fig\ire 7 >- Subjective rating probability distribution. 





Figure 10.- Satisfaction probability distribution - Twin Otter, aft seat location 




Figure 11.- Satisfaction probability distribution by aircraft type 



